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We report a series of numerical simulations showing that the critical magnetic Reynolds number 
Rmc for the nonhelical small-scale dynamo depends on the Reynolds number Re. Namely, the 
dynamo is shut down if the magnetic Prandtl number Prm = Rm/Re is less than some critical 
value Prm.c ^ 1 even for Rm for which dynamo exists at Prm > 1- We argue that, in the limit of 
Re 00, a finite Prm.c may exist. The second possibility is that Pr^.c ^ as Re oo, while Rmc 
tends to a very large constant value inaccessible at current resolutions. If there is a finite Prm,c, 
the dynamo is sustainable only if magnetic fields can exist at scales smaller than the fiow scale, i.e., 
it is always effectively a large-Pfm dynamo. If there is a finite Rmc, our results provide a lower 
bound: Rmc > 220 for Prm < 1/8. This is larger than Rm in many planets and in all liquid- metal 
experiments. 

PACS numbers: 91.25.Cw, 95.30.Qd, 47.27.Gs, 47.65. -l-a 



The simplest description of a conducting fluid is in 
terms of equations of magnetohydrodynamics (MHD): 

9tu + u-Vu = z/Au- Vp + B- VB-Hf, (1) 
StB + u-VB = B-Vu + T^AB, (2) 

where u is velocity, B is magnetic field, f is the external- 
force density, v is viscosity and ry is magnetic diffusivity. 
The pressure gradient Vp is determined by the incom- 
pressibility condition V • u = 0. We have rescaled p and 
B by p and {AirpY/"^, respectively (p is density). 

A fundamental property of Eqs. (^12) is the ability 
of the velocity and magnetic fields to exchange energy. 
In three-dimensional turbulent flows (and in many other 
chaotic flows), this can take the form of net ampliflca- 
tion of magnetic field with time, a process referred to as 
MHD dynamo. There are two kinds of dynamo. The first 
is the mean- field dynamo (growth of (B)), which usually 
requires a net flow helicity 0]. It is a large-scale effect 
that must be considered in conjunction with such system- 
specific properties as geometry, rotation, mean shear etc. 
The second kind is the small-scale dynamo: amplification 
of the magnetic energy (B^) due to random stretching of 
the field by the turbulent flow, requiring no net helicity 

IE 9 ' The stretching is opposed by the resistive diffu- 
sion, so the dynamo is only possible when the magnetic 
Reynolds number Rm = (u^)^/^£o/^ exceeds a certain 
critical value (-^o is the system scale) . 

In this Letter, we study the existence of small-scale 
dynamo in homogeneous incompressible turbulence with 
magnetic Prandtl number Prm = v/rj < 1 (i.e.. Re — 
{u^y/^£o/i' > Rm). This is an important issue because 
Prni is small in stars (Prm ^ 10^^ at base of the Sun's 
convection zone), planets (Prm^ 10~^ |5|), and in liquid- 
metal laboratory dynamos 



In three dimensions, most types of turbulence at scales 
much smaller than the system size are predominantly 
vortical and well described by Kolmogorov's dimensional 
theory In this theory, the fastest fleld stretching is 
done by the small-scale velocities. The essential physics 
of small-scale dynamo should thus be contained within 
our homogeneous, isotropic, incompressible model. 

The small-scale dynamo is most transparent in the 
limit of Prm ^ 1. Straightforward estimates show that, 
while velocity is dissipated at the viscous scale £i, ^ 
magnetic fleld can occupy smaller scales down 
to the resistive £ri ^ Vr^^^in. The dynamo is driven 
by the fastest eddies: the viscous-scale ones, which are 
spatially smooth. The growing flelds are organized in 
folds, with direction reversals at the resistive scale Iri 
and fleld lines remaining approximately straight up to 
the flow scale US O [H [11. Why this is a win- 
ning conflguration is best seen on the example of a lin- 
ear velocity fleld (locally uniform rate of strain) 0, 0| : 
the fleld aligns with the stretching direction of the flow 
but reverses along the "null" direction, so that compres- 
sion cannot lead to resistive annihilation of antiparallel 
flelds canceling the effect of stretching. For this mecha- 
nism to apply, it is essential that (i) the flow be spatially 
smooth, so fluid trajectories separate exponentially in 
time leading to exponential stretching, (ii) a scale sep- 
aration between the field scale (reversals) and the fiow 
be achievable. The large-Prm turbulent dynamo satisfies 
both conditions, as do all deterministic chaotic dynamos 
J:J. Thus, the small-scale dynamo, as it is usually un- 
derstood, is essentially the large-Prm dynamo. The often 
simulated case of Prm = 1 belongs to the same class: the 
magnetic energy is amplified at scales somewhat smaller 
than the viscous scale and the field structure is similar 
to the case of Prm > 1 (l^. When Prm < 1, the field 




FIG. 1: (a) Evolution of magnetic energy {B'^){t) in all runs. Also shown are the versions of Runs A3 and B4 that started from 
the saturated state of Run Al (at t = 75) and Run Bl (at t — 50), respectively. Runs Al, A2, and Bl have resolution 128'', 
Runs A3, B2, B3, and B4 are 256'^. (b) Growth rates vs. Prm (at fixed rf) for the same runs plus for 7 runs (128"') extending 
Series B to Prm > 1. The dynamo is again shut down at Prm ~ 100 because velocity is strongly damped by large viscosity and 
Rm drops below critical. We show Rm (divided by its value for Run Bl, Rm ~ 210) in the same plot. 



scale is resistively limited to be comparable to, or larger 
than, the viscous cutoff. The field interacts with inertial- 
range motions, which are spatially rough and cannot be 
thought of as having a locally uniform rate of strain. Is 
there still a small-scale dynamo? 

In order to address this question, we have carried out 
a series of numerical simulations. Equations were 
solved in a triply periodic box by the pseudospectral 
method. We used a random nonhelical forcing f applied 
at the box scale and white in time. The average injected 
power e = (u-f) was kept fixed. The code units are based 
on box size 1 and t — \. Defining Rm = {u^y/"^ /r]ko with 
fco — 27r the box wave number, we have found that the dy- 
namo existed for Prm > 1 provided Rm > 80 (cf. Fig.^). 
Once the existence of the dynamo at Prm = 1 for a given 
T] was ascertained, 77 was fixed. The viscosity 1^ was then 
decreased. The results are summarized in Fig. ^ In runs 
with 7] = 2 X 10^^ (series A, Rm ~ 110), starting from 
weak field, the dynamo persisted at Prm — 1/2, but was 
shut down at Prm = 1/4. In runs with 77 = 10~^ (series 
B, Rm = 220 ± 10), there was dynamo at Prm = 1/2 and 
1/4, but not at Prm = 1/8. Figure |21 shows the time- 
averaged normalized magnetic-energy spectra for series 
A and B, as well as velocity spectra multiplied by k^. 
The latter characterize the turbulent rate of strain and 
peak at the viscous scale. We see that as this scale drops 
below the resistive scale, the dynamo shuts down. Note 
that there is no initial-condition dependence: Runs A3 
and B4, which decay starting from weak field, also de- 



cay if initialized in the saturated state of their Prm = 1 
counterparts (Fig.^). 

Our main result is, thus, that Prm.c exists even as Rm 
is kept approximately fixed at a value for which small- 
scale dynamo is possible at larger Prm (Fig.^^)- In other 
words, the critical magnetic Reynolds number for growth 
Rmc increases withKe. Because of resolution constraints, 
we cannot afford a parameter scan to produce the depen- 
dence Rmc (Re). What a priori statements about this 
dependence can be made on physical grounds? 

Consider first the asymptotic case Re ^ Rm ^ 1. 
The resistive scale then lies in the inertial range, £0 ^ 
^77 ^ ^1/- As, in Kolmogorov turbulence, uijl ^ l""^/^, 
most of the stretching is done by the eddies at the re- 
sistive scale in ~ Rm~^^'^£n, where stretching is of the 
same order as diffusion [ij, |l5| . Since the inertial range 
is self-similar, the existence of the dynamo should not 
depend on the exact location of 1^ and it is the local (in 
k space) properties of the turbulent velocity that deter- 
mine its propensity to amplify magnetic energy. There- 
fore, if the dynamo fails, it does so at all Prm below some 
critical value Prm,c of order unity. The effective transi- 
tion from the "large-Prm" to the "small-Prm" regime oc- 
curs at Prm — Prm,c- In this case, Rmc/Re — > Prm.c = 
const < 1 as Re — > 00. Thus, if our results are asymp- 
totic, then the turbulent small-scale dynamo is always, 
in essence, a large-Prm one, and the folded direction- 
reversing fields are the only type of magnetic fluctuations 
that can be self-consistently generated and sustained by 
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FIG. 2: Magnetic-energy spectra normalized by {B'^)/2 and averaged over time: (a) Series A, (b) Series B. Also given are 
velocity spectra multiplied by and the reference Kolmogorov slope k^^^. The time intervals used for averaging are for 
Runs Al, A2: 5<t< 40; for Run Bl: 2.5 < t < 17.5; for Runs A3, B2, B3, B4: 10 < t < 25. 



nonhelical turbulence. As the separation between parallel 
and transverse scales of the field diminishes at Prm < 1 
(Fig. |3J), no steady fluctuation level can be maintained. 

The second possibility is that Rmc asymptotes to some 
constant value for Re above those we are able to resolve: 
Kmc — > const > 220 and Prm c ^ as Re — > oo. Our 
results do not rule out this outcome, whereby asymptoti- 
cally in Rm and Re, the dynamo persists at low Prm, but 
very large Rm is needed to achieve it in practice (numer- 
ically or experimentally). In stellar convective zones, Rm 
is, indeed, very large (10^. ..10^ for the Sun). On the other 
hand, in planets and in laboratory dynamos, Rm ^ 10^, 
which is comparable to Rm in our simulations. 

Note that the arguments above assume scale invari- 
ance of the inertial range, i.e., neglect the effects of in- 
termittency. An intermittent velocity field will exhibit 
large coherent fiuctuations of the rate of strain, which 
might be locally effective in stretching the magnetic field 
in a way similar to the large-Prm dynamo . Whether 
these fiuctuations can provide enough stretching on the 
average to make a workable dynamo cannot be settled 
qualitatively. Note that an intermittent growth by rare 
strong bursts is evident in Pr^ < 1 runs where the dy- 
namo is suppressed but not shut down (most vividly in 
Run A2, see Fig.[lji). 

No theory of the dynamo shutdown at low Prm ex- 
ists at present. Invoking turbulent diffusion (mixing) of 
magnetic fields by the subresistive-scale motions as the 
suppression mechanism makes heuristic sense, but does 
not provide an unambiguous verdict on the existence of 
the dynamo. Indeed, when £o ^ £ri ^ i^, the dominant 
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FIG. 3: Characteristic scales for Series B (cf. Fig.^t): fcrms = 
((j VB|^)/(_B^)) ^''^ is (roughly) the inverse direction-reversal 
scale, = ((|B • \7B\^) / {B*))^^'^ is the inverse fold length. 



and kx = {{\Vu\'^) / (u^)) ' is the inverse Taylor microscale 
(times \/5) (see R,efs. ll^ll3l for discussion of these quantities). 



contributions to both stretching and mixing arc from 
the resistive scale £ri, and the outcome of their compe- 
tition is impossible to predict on a heuristic basis. Many 
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aspects of small-scale dynamo have received successful 
theoretical treatment in the framework of the Kazant- 
sev model which assumes a Gaussian white- in-time 
velocity field, x)M-J(t', x')) = S{t - t')K'^{yL - x'). 

The correlator can be expanded, K*-'(y) — k*-'(0) ~ ~J/", 
when y ^ £b, the magnetic-field scale. When Pr^ 3> 1, 
<C iu, so a — 2, corresponding to the spatially 
smooth viscous-scale eddies. On the other hand, when 
Pr,,! 1, £b ^ ^ and magnetic field interacts 
with rough inertial-range velocities, which must be mod- 
eled by a < 2. The Kazantsev velocity is not a dynamo if 
it is too spatially rough, viz., when a < 1 | la llTl llsL Il9| . 
If we set aside fundamental objections to the (5-correlated 
model and try to compare it to real turbulence, we still 
face the difficulty of interpreting the 6 function. If we 
write equal-time velocity correlators by replacing the S 
function by inverse correlation time 1/tc, then the rela- 
tion between a and the spectral exponent of the turbu- 
lence depends on how we choose Tc. The usual choice for 
Kolmogorov turbulence is ~ y'^^^, the eddy-turnover 
time. Then a ~ 4/3 > 1 and there is dynamo [l5l|[2fil|. 
Note that Rm ^ in this case is typically much larger than 
for a = 2 [l^, lis I2OI ■ Although specific values of Rmc 
calculated from the Kazantsev model cannot be consid- 
ered as quantitative predictions for real turbulence, they 
appear to point to the second possibility mentioned above 
(finite but unresolvably large Rmc). We emphasize that 
all these results depend on the heuristic choice of Tc (e.g., 
if Tc const, a = 2/3 and there is no dynamo) and on 
the universality of the physically nonobvious condition 
a > 1. It is fair to observe that our simulations are still 
too viscous to have a well-developed Kolmogorov scaling 
(Fig. 121). Thus, if the existence of the dynamo depends 
on the exact inertial-range scaling of the velocity field 
and/or only manifests itself at very large Rm, neither 
the Kazantsev theory nor simulations at current resolu- 
tions can lay claim to a definitive answer. Obviously, the 
Kazantsev theory also cannot capture any role the inter- 
mittency of the velocity field might play and, more gener- 
ally speaking, it is doubtful that a (5-correlated Gaussian 
flow is a suitable model of the inertial-range turbulence. 

While, as far as we know, ours is the first systematic 
study of the small-scale dynamo suppression in homoge- 
neous isotropic MHD turbulence with low Pr^, indica- 
tions of this effect have been reported in the literature in 
two previous instances. Dynamo suppression at low Pri„ 
was seen by Christensen et al. ^1] in their simulations of 
convection in a rotating spherical shell and by Cattaneo 
in simulations of Boussinesq convection. These cases 
of failed low-Prm dynamos in inhomogeneous convection- 
driven turbulence are likely to be related to the same 
universal mechanism that made the dynamo inefficient in 
our simulations. Our key conclusion is that the dynamo 
suppression is a generic effect unrelated to the particular 
type of driving or other large-scale features of the system. 

The mean-field dynano, which, in contrast, does de- 



end on large-scale features such as helicity and rotation 
1| , may then be the only type of self-sustained field gen- 
eration for low-Prm systems. If a mean field is present, 
it gives rise to a source term (B) ■ Vu in the induc- 
tion equation Q and thus induces small-scale magnetic 
fluctuations. They have a fc~^^/^ spectrum at k ^ k 
[1^ I15I l23j . which has been seen in the laboratory 
and in large-eddy simulations jl^l . 
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